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Motivated by the experimental accessibility of rare B decays in the ongo- 
ing and planned experiments, we propose to undertake a model-independent 
analysis of the inclusive decay rates and distributions in the processes 
B — ► Xgj and B — > X s £ + £~ (B = 5 ± or B®). We show how measure- 
ments of the decay rates and distributions in these processes would allow 
us to extract the magnitude and sign of the dominant Wilson coefficients of 
■ the magnetic moment operator m^s^o iW bnF ilv and the four-fermion opera- 

tors (sl7^l)(^7 / ^) and (sl7 A4 6l)(£7 a '7 5 £). Non-standard-model effects could 
thus manifest themselves at low energy in rare B decays through the Wil- 
son coefficient having values distinctly different from their standard-model 
counterparts. We illustrate this possibility using the examples of the two- 
doublet Higgs models and the minimal supersymmetric models. The dilepton 
invariant mass spectrum and the forward-backward asymmetry of £ + in the 
centre-of-mass system of the dilepton pair in the decay B — > X s £ + £~ are 
also worked out for the standard model and some representative solutions 
for the other two models. 
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1 Introduction 



The measurement of the decay mode B — > K*^y by the CLEO collaboration 
[|], having a branching ratio £>(£? — > K*7) = (4.5 ± 1.0 ± 0.9) x 10~ 5 , has 
put the physics of the electromagnetic penguins on an experimental foot- 
ing. This measurement and the experimental upper bound on the inclusive 
decay BiB — > X s ^) < 5.4 x 10~ 4 at 90% C.L. g have been analysed in 
the context of the standard model (SM) || and in extensions of it such as 
the two-Higgs-doublet models (2HDM) || [5|, the minimal supersymmet- 
ric models (MSSM) |7|, §, |§, and a number of other more exotic variations 
[ITD, [TTj| . Very recently, the CLEO collaboration has published a first mea- 



surement of the inclusive decay rate, based on the measurement of the photon 



energy spectrum in the decay B — > X s ^f |12J . The branching ratio 



B(B -»■ X sl ) = (2.32 ± 0.51 ± 0.32 ± 0.20) x 10" 4 (1) 

puts more restrictive bounds on the non-SM parameters. In the SM context, 
the short-distance contributions in these transitions are dominated by the top 
quark and hence they provide valuable information about its mass and the 
Cabibbo-Kobayashi-Maskawa (CKM) weak mixing matrix elements V ts V t b 
[ fn| . Specifically, the recent CLEO measurement of E{B — > X s ^) yields at 
present the following bounds M: 



0.62 < 



f.s 



V, 



cb 



< 1.1 • (2) 



Alternatively, the ratio V ts /V c b can be determined from unitarity and one 
may use | Vj s |/| V c ft| ~ 1 to obtain from the CLEO measurements bounds on 
the Wilson coefficient Cj{rrib) of the effective magnetic moment operator. 
Using B(B -> X s 7) = (2.32 ±0.67) x 10~ 4 , obtained by adding the statistical 
and the systematic errors in quadrature, one obtains 

0.22 < \C 7 (m b )\ < 0.30. (3) 

This bound is subject to the residual next-to-leading order corrections. The 
inclusive branching ratio B(B — > X^) = (2.32 ± 0.67) x 10~ 4 is consistent 
with the exclusive branching ratio B(B — > K* + 7) = (4.5 ± 1.0 ±0.9) x 10~ 5 , 
with R(K*/X S ) = T(B -> K* + i)/T(B -> X s + 7) calculated to be 0.1 < 
R(K*/X S ) < 0.2 in most theoretical models of recent vintage 0, fl4]|-|^ . 
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The bound on | CV(m&)| given in ([3[) can be used to constrain the non-SM 
contribution to the decay rate BiB — > -Xyy). The inclusive branching ratio 
in the SM, including the leading corrections in the anomalous dimension ma- 
trix and 0(a s ) virtual and bremstrahlung corrections, has been estimated 
to be B(B -> X sl )= (2.8 ± 1.0) x 10~ 4 @. With m t = 174 ± 16 GeV de- 
termined from the CDF data |22|, a good part of the uncertainty in the 
theoretical estimates of BiB — > X s ^) is due to the scale dependence of 
the perturbative QCD framework and the QCD scale parameter entering 
ag itself. Recently, parts of the next-to-leading-order corrections have been 
included in the anomalous dimension matrix and the resulting decay rate 



T(B — > X s + 7) has been recalculated [^TJ. The inclusion of these terms 



reduces the scale-dependence of the effective coefficient CV(/x) as the scale fi 
is varied in the range m^jl < p < Im^. The partial next-to-leading-order, 
however, not unexpectedly shows a regularization scheme-dependence and 
its value is also somewhat diminished compared to the leading-log estimates. 
This reduces the branching ratio, and the (partial) next-to-leading-order es- 



timates give B(B -> X s 7)= (1.9 ± 0.2 ± 0.5) x 10" 4 J21|. The complete 
next-to-leading-order result is not yet available and hence the estimates of 
the SM for BiB — > X s '-f) are not completely quantitative, although the in- 
clusive rate (P and the SM estimates || ^l|] are in good agreement. The 
sensitivity to new physics in B — > X s ^ is therefore somewhat entangled with 
the QCD corrections. Despite this, the CLEO data on rare B decays have 
provided valuable constraints on the parameters of models that are popular 
candidates for the extensions of the SM, and in some cases these constraints 
are very competitive compared to the ones following from direct and indirect 



searches [12 



The determination of \C-j(jrib)\ from the inclusive branching ratio B(B — > 
Xgj) is a prototype of the kind of analysis that we would like to propose here 
to be carried out for the rare B decays in general and for the semileptonic 
decays B — > X s £ + £~ , in particular. In the standard model, a measurement 
of the radiative rare B decay B — > X^ mentioned above and related ones, 
such as B — y X s £ + £r, B — > Xa the corresponding exclusive decays 

B — > (K, K* , 7r, p, ...)£ + £~ , would lead to a determination of m t and the 
CKM matrix element product V ts V t b and V t dV t b- The purely leptonic decays 
B° s £+£- and B° d -> £ + £~ and the modes B° s 77, B Q d 77, while hav- 
ing very different final states and branching ratios - and hence representing 
experimentally very different propositions - have more or less the same infor- 
mation content in the SM as the decays mentioned above (see, for example, 
and references cited therein). 

The aim of this paper is to undertake first steps towards a model-independent 
analysis of the FCNC electroweak rare B decays. Although the method of 
analysing the data and the relevant theoretical framework being presented 
here are developed for the inclusive decays B — > X s ^ and B — > X s + £ + £~ , 
much of the general considerations being discussed apply also to the corre- 
sponding exclusive decays such as B — > K*j, B — > K* £ + £~ and B — > 



K£ + £~ . Of course, the extraction of the short-distance physics in terms of 
the Wilson coefficients of the dominant operators from the data on exclu- 
sive decays would require the knowledge of the relevant form factors. This 
may compromise the precision on the short-distance part of the amplitudes 
in question at present, though advances in computational methods for QCD 
may allow quantitative conclusions to be drawn from exclusive decays also. 

Our analysis is based on the renormalization group (RG)-improved treat- 
ment of the effective Hamiltonian relevant for B decays under consideration, 
obtained by integrating out the top quark and the other heavy degrees of 
freedom. The resulting Hamiltonian in the SM can be written as: 

AC 10 

H eff (b - sX) = -r^GMOiM . (4) 
V2 i= i 

Here X stands for qq, 7, gluon and £ + £~ and X t = V t *V t b] the operator basis 
contains dimension-5 and -6 operators and is given in the appendix. 

The problem of carrying out a model-independent analysis in FCNC pro- 
cesses is evident from the SM expression for Ti. e ff given above. In the most 
general case, which would also include Left-Right-symmetric (LR) models, 
the operator basis consists of 20 operators, having 20 independent coeffi- 
cients Ci. However, as argued later, we shall limit the present analysis to 
left-handed fields only, in which case there are still 10 operators to be consid- 
ered. The effective Hamiltonian then still involves 10 independent coefficients 
Ci and determining them experimentally is a monumental task. It is obvious 
that some theoretical assumptions have to be used to focus the attention 
on the more interesting Wilson coefficients. In our opinion, the promising 
coefficients from the stand-point of non-SM physics searches are the ones in 
which the electroweak loop effects (penguins and boxes in the diagrammatic 
language) play the dominant role. 

The analysis carried out here is restricted to the models in which the effec- 
tive Hamiltonian is of the form given in (|4]). We show later how measurements 
can be analysed to determine the inadequacy of this operator basis, should 
that happen. Although by no means completely general, the operator basis 
contained in W e // encompasses most models of current theoretical interest, 
which include, apart from the SM, the 2HDM (both types I and II) and the 
MSSM. The effective Hamiltonian in the four-generation models, as well as 
in models with anomalous trilinear gauge and fermion couplings, can also be 
written in the above form. 

The coefficients Ci(Myy),i = 1, . . . , 6, determine the non-leptonic B- 
decay rates and the .B-hadron lifetimes. Since data on B decays and the 
results obtained in the SM from (|]), including QCD effects, are in good 
agreement with each other (±20%), there is not much room left for the first 
six coeflficients involving the four-quark operators to deviate from their SM 
values. So we fix Ci(M w ),i = 1,...,6, to their SM values. Beyond the 
leading- log approximation, these receive corrections of order a s (M w ) at the 
large scale /1 = M w ; they should be incorporated in theoretical estimates 



consistently with other higher-order effects. The coefficients of our interest 
are CV(/i), Cio(fi), since they are generated at scale \i = Mw by elec- 
troweak loops (penguins and boxes). They govern the physics of the rare B 
decays b -> s + g, B -> X s7 , B -> X s and 5° -> We shall con- 

centrate on them and show how to extract these coefficients (both their signs 
and magnitudes) from data on radiative and (semi)leptonic rare B decays. 
These can be compared with the SM and extensions of it to search for new 
physics. 

The experimental quantities we consider in this paper are the following: 

• Inclusive radiative rare decay branching ratio B(B — > -XVy); 

• Invariant dilepton mass distributions dB(B — > X s £ + £~)/ds; 

• Forward-backward (FB) charge asymmetry A(s) in B — > X s £ + £~ . 
The kinematic variables are defined as: 

u = (p b -Pif - (Pb -P2) 2 , 
s = ( pi +p 2 ) 2 ; 
s 

S 2 ' 

mi 



w(s) = \J(s - (m b + m s ) 2 )(s - (m b - m s ) 2 ). (5) 

where p b , P\ and p 2 denote, respectively, the momenta of the b quark (= B 
hadron), £ + and £~ . The FB asymmetry A(s) is defined with respect to 
the angular variable z = cos 6 = u/w(s), where 9 is the angle of the £ + 
with respect to the 6-quark direction in the centre-of-mass system of the 
dilepton pair. It is obtained by integrating the doubly differential distribution 
d 2 B/(dzds) H: 
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The rationale for concentrating on these measurements is the following. 
We remark that the decay rate B{B — > X s7 ) puts a bound on the absolute 
value of the coefficient Cj(fi). However, the radiative B decay rate by itself 
is not able to distinguish between the solutions C 7 (fi) > (holding in the 
SM) and the solutions Cr(/j) < 0, which, for example, are also allowed in the 
MSSM as one scans over the allowed parameter space. We recall that the 
invariant dilepton mass distribution and the forward-backward asymmetry 
in B — » X s £ + £~ are sensitive to the sign and magnitude of C 7 (/i) [fZ|, gH ]. 
This is easy to see in the approximation of neglecting the s-quark mass, in 
which limit the differential branching fraction for B — > X s ^ can be written 
as: 

dB{B " d f ^ = K(1-S) 2 {(|^( / ,)| 2 + |C 10 (^)| 2 )(1 + 2S) 

+|C T Gu)| 2 ^(2 + s) + 12 Re (C 7 (/i)C 9 e// (/i))}, (7) 
s 



where K is a constant and Cg (//) is defined as 

C e 9 ff (fi) = C 9 (fi) + Y(^s) , (8) 

where Y(fi, s) is a function involving the coefficients Cj(/i); i = 1, . . . , 6, and 
it depends on the kinematic variable s through the one-loop matrix element 
of the four-quark operators. Long-distance effects would also be contained 
in Y(fi, s); however, we shall not consider these here. 

Consistent with our assumptions, we shall take the function Y(fi) from 
the SM and admit non-SM contributions only in CV(/i), Cg(/i). Cg(jj), and 
Ciq. From the above expression it is obvious that the dilepton invariant-mass 
distribution is sensitive to the real part of the product CV(/i)Cg (/i) in sign 
and magnitude. However, this distribution by itself cannot determine the 
sign of Cg^ and C\q. As the FB asymmetry and the dilepton rate are in 
general independent quantities, they provide independent constraints on the 
Wilson coefficients, and we argue that the FB asymmetry A(s) can be used 
to resolve these ambiguities. In the limit m s = 0, the FB asymmetry has a 
simple form also: 

A(s) = ~K(1 - S)*CM{<%"W + ^MY (9) 

The s-dependence in A(s) now allows us to disentangle the dependence on 
CV(/i) and Cg (/i). Moreover, the FB asymmetry A(s) is directly propor- 
tional to the coefficient Cio(/-0 an d hence the asymmetry is very effective in 
constraining Ci (/i). 

The first measurements will be done for partially integrated dilepton mass 
distribution and asymmetry. We discuss how to extract the Wilson coeffi- 
cients from these measurements. A given integrated branching ratio B(As) 
over a given range in s, As, and a given value of CV(/x) can be analysed in 
terms of the contour plots drawn in the Wilson-coefficients plane Cg(fj,) and 
C w (/j,). It is obvious from the quadratic equation (^) that, for a given C 7 (fi), 
there is at most a four- fold ambiguity in the signs of C 9 (/x) and C w (fi). In 
addition one may make use of the integrated FB asymmetry A(As). For 
a fixed value of A(As) one finds hyperbolae in the Cq-C w plane, for given 
sign and magnitude of CV(/z). The intersection of the contours obtained from 
B(B — > X s £ + £~; As) and A(As) then determines the signs and magni- 
tudes of all three Wilson coefficients C-i(fi),Cg(fi) and Cio(/x). We plot, in 
the Cg(ju)-Cio(ju) plane, the contours that are allowed from the solutions of 
the two algebraic equations for 13 (As) and A(As) for the two values of C 7 , 
Cj = ±0.3, and two invariant-mass intervals below the J/ip and above the ip' 
mass (see figs. pointing out the SM solution and the general solutions in 
non-SM cases. It is, however, conceivable that there is no consistent solution 
of the two algebraic equations for 13 (As) and ^4(As) in B — > X s f + £~ and 
the branching ratio B — > X s 7 in terms of the coefficients Cj(/i), Cg(n) and 
C w (fi), which would then indicate that the operator basis chosen in describ- 
ing the Ti-eff given above is incomplete. The LR-symmetric models, which 



we have dropped from our discussion, is a case in point. Of course, given 
enough statistics the two distributions could be used to determine the coeffi- 
cients from the data fits and compared with the values in the various models 
directly. We illustrate these analysis techniques for some assumed values of 
the branching ratios and asymmetry and representative values in the allowed 
parameter space in various models in terms of both the contour plots in the 
Wilson coefficient space and the distributions themselves. 

This paper is organized as follows. In section 2, we briefly review the 
kinematics of the decays B — > X s ^ and B — > X s £ + £~ , define the amplitudes 
and the various distributions of interest. In section 3, we perform a numerical 
study of the differential and integrated branching ratio in B — > X s £ + £r . To 
illustrate the method, we have taken a representative value for the effective 
Wilson coefficient, ^(mj) = ±0.3, which lies within the presently allowed 
range for this quantity from data on B — > X^. We give the partially in- 
tegrated branching ratio B(As) and the asymmetry A(As), making the de- 
pendence on the coefficients Cg and C\q explicit. The resulting constraints 
on these coefficients from the "low-dilepton mass" and "high-dilepton mass" 
regions are then displayed as contour plots in the Cg-Cw plane for some 
illustrative values of the branching ratios and FB asymmetries. The invari- 
ant dilepton mass spectrum in B — > X s £ + £~ and the FB asymmetry are 
also shown here for some representative values of the Wilson coefficients. In 
section 4 we discuss the predictions of specific models, namely the 2HDM 
(type I and type II) and the MSSM. These models allow a calculation of the 
coefficients CV(/i), C 9 (/i) and Ci (/i) as a function of the model parameters. 
For the case of the 2HDM we show CV(/i), C 9 (/x) and Ci (/i) as a function of 
the charged Higgs mass for representative values of the ratio of the two 
vacuum expectation values v^jvy. As already noted [^|, the type-I 2HDMs 
admit negative CV(/u) solutions, as opposed to the positive Cj(fi) solution 
obtained in the SM. However, for the type-I models the negative C 7 (/i) solu- 
tions are excluded by present data, as will be shown in section 4. To analyse 
the MSSM we vary the parameters of the model over the experimentally al- 
lowed values, and the resulting region in the Cg(/i)-Cio(/u) plane is shown. 
In doing this, we have imposed the constraints on the coefficient CV(/i) as 
explained earlier. The restrictions on the supersymmetric (SUSY) particle 
masses from present searches and the anticipated reach of experiments at 
the Tevatron and LEP-II are also imposed. The most interesting part of 
this exercise is that these constraints do allow both the negative and positive 
CV(/x) solutions, which are consistent with the data on B — > X s ^ but admit 
values of the coefficients Cg(/i) and C*io(/i) sufficiently different from those of 
the SM, thereby yielding very different differential distributions in the decay 
B — > X s £ + £~ . These distributions are shown for some representative values 
of the Wilson coefficients. Section 5 contains a summary. To fix our notation 
we collect the relevant formulae for the effective Hamiltonian, including the 
RG-improved Wilson coefficients Cj(/i) in an appendix. 



2 Estimates of B(B — > X 8 j), Invariant Dilep- 
ton Mass Distribution and Forward— Backward 
Asymmetry in B —> X s £ + £~ 

In the following we shall consider the inclusive decays B — > and B — > 
X s £ + £~, where £ is either electron or muon. It has been shown that inclusive 
B decays may be treated in a 1/m^ expansion, the leading term of which 
is the free quark decay. The next-to-laeding effects are of second order in 
i.e. 0(m 2 s /ml) 127], ^8], E9]. This is true for total rates, for partially 
integrated rates, and also for decay distributions, as long as one is not too 
close to the kinematic endpoint in which the energy release in the hadronic 
subprocess becomes small. 

We shall also not include QCD corrections to the free-quark decay distri- 
butions, aside from the leading logarithms, which are induced by the renor- 
malization group running. It is known that the measured inclusive lepton 
energy distributions of the charged- current- induced semileptonic transitions 
are well approximated by the partonic distributions b — > (c, u)£~v^ including 
QCD corrections |TT], [32J. We expect this for the flavour-changing neutral 
current semileptonic transitions as well. 

The second proviso is that we shall concentrate on the short-distance 
contributions to the decays B — ► X s ^ and B —>■ X s £ + £~ . This again is not a 
drastic oversimplification, as it is known that the long-distance contributions 



to the decay B — ► X s "f are small [33, B3, 35 1 and those in the decays B 



X s £ + £~ are dominantly present at and near the J / ip and ip 1 poles, extending 
also to the region between them. These contributions can be modelled using 
data on the decays B — > (J/ip,ip')X s — > (£ + £~)X s , which already exist and 
which are expected to become quite precise in the future. We shall thus 
concentrate on the distributions in the regions away from the resonances and 
consider the two regions 

"Low dilepton-mass" : 4m 2 t < s < mju — 5 , 

"High dilepton- mass": m^, + S < s < s max , (10) 

where 5 is a cut-off that can be matched with the experimental cuts used 
in the analysis. In these two regions we expect only small long-distance ef- 
fects. Due to limitations in rates and dilepton trigger requirements, we antic- 
ipate that the "Low dilepton-mass region" is accessible to e + e~ experiments 
(CLEO and B factories) and the "High dilepton-mass region" typically to B 
experiments which will be carried out with hadron beams (CDF, HERA-B, 
LHC). Given high enough luminosity, this region can also be probed at the 
B factories. 



2.1 Decay rate for B —> X s "f 



The procedure for the computation of the Wilson coefficients in the SM and 
extensions of it is by now standard. They are obtained at a large scale 
/i 2 = Myy (in general at /i 2 = m 2 ^, where m Y stands generically for the mass 
of the heavy degrees of freedom) by integrating out degrees of freedom heav- 
ier than jj, P^| , and are sensitive to the presence of new physics. Then the 



renormalization group (RG) equations are used to scale down these coeffi- 
cients to the scale that is typical for i?-hadron decays, namely fi = 0(m b ). 
In this way large logarithms of the form a s In are shifted from the ma- 
trix element into the Wilson coefficients. The RG evolution of the Wilson 
coefficients Cj(/x) involves the (10 x 10) anomalous dimension matrix, which 



has now been calculated at the one- loop level [37], |38|. This matrix and the 
solutions of the RG equations in terms of the effective coefficients Cj(/i) are 
given in the appendix. 

Starting from the effective Hamiltonian as given in the appendix, one finds 
that only one operator, namely O7 contributes at tree level. The contribution 
of the QCD bremsstrahlung process b — > 57 + g and the virtual corrections 
to b —>■ sj have been calculated in 0(aa s ) in ref. f3"0| . Expressed in terms of 
the inclusive semileptonic branching ratio BR(B — > Xli/i), one can express 
the branching ratio BiB — > Xg'y) as: 



B(B^X sl ) = 6- / u 1 V] xB sl , (11) 

7r \V cb \ z f{m c /m b )[l - [2a s )/{37r)h{m c /m b )\ 

where 

f{r) = 1 - 8r 2 + 8r 6 - r 8 - 24r 4 ln(r) 

is the phase-space function for Y{b — > c + ivi) and B s i = 10.5% is the semilep- 
tonic branching fraction. The function h(r) accounts for QCD corrections 
to the semileptonic decay and can be found, for example, in ref. [3~I]. It 



is a slowly varying function of r and, for a typical quark-mass ratio of 
r = 0.35 ± 0.05, it has the value h(r) = 2.37 =F 0.13. The contributions 
from the decays b — > u I vi have been neglected in the denominator in (|TT|) 
since they are numerically inessential because \y ub \ IK&I- 

The inclusive decay width for B — > is dominantly contributed by 
the magnetic-moment operator Or, hence the rationale of factoring out its 
coefficient in the expression for B(B — > Xsj) in Eq. (|TT|). Including 0(a s ) 
corrections brings to the fore other operators with their specific Wilson co- 
efficients. The effect of these additional terms can be expressed in terms of 
the function K(n), which lumps together the effects of bremsstrahlung cor- 
rections in the inclusive decay rate. The function K([i) has been computed 
in the SM in [|J taking into account the dominant corrections from C2 and 
C$ (the coefficients of other operators are considerably smaller). Typically, 
0.79 < K(ji) < 0.86 for m b /2 < jj, < 2m b . Since C 8 (/i) receives possible 
non-SM contributions in the same way as CV(/i), one should compute the 
K-factor and incorporate these corrections into a full next-to-leading-order 



analysis. However, we are not working beyond leading order in this paper 
and so we ignore the K factor. Using the CLEO inclusive measurement 



B{B -> X sl ) = (2.31 ± 0.67) x 1(T 4 , (12) 

we get the following bound on CV(/^): 

0.22 < |C 7 (//)| < 0.30. (13) 

Using, however, the 90%-confidence-level range from the CLEO measurement 
B(B -> X s 7) = (2.31 ±1.1) x 10" 4 and the theoretical calculation for B(B -> 
X s 7) from |J we obtain 

0.19 < |C 7 (//)| < 0.32. (14) 

This range, in our opinion, adequately reflects the present uncertainties. 
When not stated otherwise, we shall fix \Cj(fi)\ = 0.3, which is in com- 
fortable agreement with the CLEO data. 

2.2 Decay distributions in B — > X s £ + £~ 

Using 7i e ff given in (|4|), one obtains for the dilepton invariant mass distri- 
bution 



2 \ 2 



dB - 

ds sl A^\V cb \ f(m c /m b ) W[S) 



(\C 9 + Y(s)\ 2 + C 2 w ) ax(s,m s )(lB) 



4 

+-C 2 a 2 (s,m s ) + 12a 3 {s,m s )C 7 (C 9 + Re Y(s)) 
s 

where the auxiliary functions are defined as follows: 



ai(s,m s ) = -2s 2 + s(l + m 2 ) + (1 - m 2 ) 2 (16) 

a 2 (s,m s ) = -(l+m2)s 2 -(l + 14m2 + m 4 )s + 2(l + m^)(l-m2)fl7) 

a 3 (s,m s ) = (l-m 2 s ) 2 -(l + m 2 s )s (18) 

Y(s) = g (m c /m b J)(3C l + C 2 + 3C 3 + C A + 3C 5 + C 6 ) (19) 

-l g (l,s)(AC 3 + AC A + 3C 5 + C 6 ) 
~g(0,s){C 3 + 3C 4 )+AC 9 

and g(z, s) is the one-loop function given in the appendix. Furthermore, the 
constant ACg depends on the scheme, which is chosen in the evaluation of 
the one- loop matrix elements of the operators 0\ ■ ■ ■ (D§. It contains also 
a large logarithm ln(M^/mj;), which is not due to QCD effets but rather 



comes from the one- loop matrix elements of 0\ ■ ■ ■ 0%. In the MS scheme 
one obtains |25 



ACq 




'a s (M) 



w i 



-11/23' 



a s (m b ) 



]7 



'a s (M } 



w ) 



-29/23' 



a s (m b ) 



which generates the correct logarithms in the limit a s — > 0. 

The corresponding differential asymmetry as defined in (|6]) is 



(20) 



A(s) 



3a 2 



In 2 f(m c /m b ) 



w 2 {s)C 10 s(C 9 + ReY(s))+4C 7 (l + m 2 pl) 



In the subsequent section, we shall be evaluating the partial branching 
ratio B(As) and partial FB asymmetry A(As), where As defines an interval 
in the dilepton invariant mass. All numerical calculations are done with a 
non-zero value for the s-quark mass, m s = 500 MeV. The asymmetry in 
the dilepton angular distribution in the SM can be qualitatively understood 
as follows. The decays B — > X s £ + £~ occur through 7 , Z and W + W~ 
exchange diagrams. For small m t (m t /M\y <C 1) the photon contribution 
dominates and the vector-like interactions to the leptonic current remain 
substantial; consequently, the asymmetry is small. However, for mt/Mw > 2, 
as suggested by the CDF value for m t , the contribution from the Z-exchange 
diagrams becomes important and the coefficient of the left-handed leptonic 
current grows as m 2 , leading to a large asymmetry. 



3 Analysis of the Decays B —> X s ^f and B —> 

In this section we shall discuss how the Wilson coefficients appearing in the 
effective Hamiltonian may be extracted from the experimental information. 
We shall assume that all the matrix elements are normalized at the scale 
fi ~ rrih, the mass of the b quark and hence the decay distributions are given 
in terms of the Wilson coefficients at the scale m&. The SM makes specific 
predictions for these coefficients (modulo perturbative QCD uncertainties), 
but if there is physics beyond the SM, these coefficients will in general be 
modified. 

We will somewhat elaborate on this point. A specific model provides the 
set of Wilson coefficients at high scales, which we shall choose to be the scale 
of the weak bosons \x = Mw- Furthermore, we shall integrate out heavy 
degrees of freedom at the same scale \i = Mw', this procedure introduces 
an uncertainty due to the difference in the masses of the heavy degrees of 
freedom, as for example arising from m t ~ 2M\y. However, since the QCD 
coupling constant is small at these very high scales and does not appreciably 
change between these thresholds, it is a reasonably accurate approximation to 



neglect QCD corrections for scales above \x = Mw- Starting from this scale, 
the Wilson coefficients are obtained from the solution of the renormalization 
group equations at the scale /i ~ m b , where we use the one-loop result for 
the anomalous dimensions and the beta function (see appendix). 

In order to determine the sign of Cj and the other two coefficients C$ 
and Cio, one has to study the decay distributions and rates in B — > X s £ + £~ , 
where £ is either electron or muon. As already discussed, these decays are 
sensitive to the sign of CV, and to Cg and C\q. The first experimental in- 
formation available in the decay B — > X s £ + £~ will be a measurement of the 
branching fraction in a certain kinematic region of the invariant mass s of 
the lepton pair. In order to minimize long-distance effects we shall consider 
the kinematic regime for s below the J/ip mass (low invariant mass) and for 
s above the mass of the ip' (high invariant mass). Integrating (20) over these 
regions for the invariant mass one finds^ 



B(As) = A(As) (C 2 + C? J + B(As)C 9 + C(As), (22) 

where A, B and C are fixed in terms of the Wilson coefficients C\ • • • Cq and 
Cj. We derive from (|HS|): 

a 2 1 r 

A{As) = B sl — £/ _ j- N / d£u(s)a x {s,m a ) (23) 
a 2 I 



As 



'47T 2 f(m c /m b ) 
a 2 1 r 

B(As) = B d —— -— / dsu(s) 2ai(s,m a Re Y(s) + 12C 7 a 3 s,m^24 



As 



As 



4 

H — |C 7 | 2 a!2(s,?fi s ) + 12a 3 (s,m s )C 7 Re Y(s 
s 



a x {s,m s ) {( Re Y{s)) 2 + ( Im Y{s)) 2 } 

(25) 



where the auxiliary functions a,, i = 1,2, 3, are as given above. 

In our analysis we keep the values for C\ ■ • ■ Cq and the modulus of CV 
fixed and hence A(As), B(As) and C(As) may be calculated for the two 
invariant-mass ranges of interest. For the numerical analysis we use m b = 4.7 
GeV, m c = 1.5 GeV, m s = 0.5 GeV. The values for the Wilson coefficients 
Ci, . . . , Cq are given in the appendix. The resulting coefficients A, B, and C 
are listed in Table [TJ. 

Inspection of ( p5|) shows that the integrand for C(As) behaves as 1/s for 
small values of s, leading to a logarithmic dependence of C(As) on the lepton 
mass for the case of the low invariant mass region. In fact, this is the only 

3 In performing the integrations over As we have set the resolution parameter S to zero, 
since we do not consider any long-distance contribution. The long-distance contribution 
peaks strongly at the J/ tp and ip 1 and 5 has to be several times the width of these resonances 
in order to avoid large long-distance effects. However, calculating only the short-distance 
part one may safely neglect S, since the short-distance contribution is flat in this region. 



As 


c 7 


A(As)/10" 8 


fi(As)/10" 8 


C(As)/10- 8 

£ = e 


C(As)/10" 8 
£ = H 


Amj < s < m 2 ^ 
4m 2 <s< m 2 ^ 
m\, < s < (1 - m 2 ) 
m\i < s < (1 - m 2 s ) 


+0.3 
-0.3 
+0.3 
-0.3 


2.86 
2.86 
0.224 
0.224 


-5.76 
-20.8 
-0.715 
-1.34 


84.1 
124 
0.654 
2.32 


76.6 
116 
0.654 
2.32 



Table 1: Values for the coefficients A(As), B(As) and C(As) for the decay 
B -> X s £+£-. 



point where the masses of the leptons enter our analysis, and from this one 
may obtain the corresponding coefficients for £ — lm: 

C{4ml < s < m 2 ^) = C(4m 2 e < s < m 2 w )-8|C 7 | 2 (l+m^)(l-m2) 3 In 

'(26) 

Of course one may apply (|26|) also to obtain C for any other lower cut sq on 




the lepton invariant mass, as long as so ^j/^,- 

For a measured branching fraction B(As), one can solve the above equa- 
tion for B(As), obtaining concentric circles in the Cg-CiQ plane, with their 
centre lying at C* = B(As)/(2A(As)) and C{ Q = 0. The radius R of these 
circles is proportional to 



R=^B(As)-B min (As), (27) 
where the minimum branching fraction 

B mm {As) = C(As) - (28) 

is determined mainly by the present data on B — > X s r y 1 i.e. by \C-j\. For the 
cases of interest one obtains, with the help of Table 0: 

„ ,, 2 2 N [ 8.1 x 10- 7 for C 7 = 0.3 /on , 

B mm {4m e < s < mj,+) = | g 6x 1Q _ 7 for Ci = _ Q 3> (29) 

2 n 2^ / 8.5 x 10- 10 for C 7 = 0.3 

tJmin{ m ^' < s < [1 rn s )) — i 3 x 10~ 9 for C 7 = — 3 

Note that B(As) is an even function of Cio, so one is not able to fix the 
sign of Cio from a measurement of B(As) alone. 

To further pin down the Wilson coefficients, one could perform a mea- 
surement of the forward-backward asymmetry A, which has been defined 
above. The asymmetry is an odd function of Cio, an d for a fixed value of the 
total branching ratio in this kinematic region one obtains, from integrating 



As 


c 7 


«(As)/10- y 


[3(As)/10- 9 


4mj < s < m 2 ^ 
4mj < s < m 2 j^ 
< s < (1 — m 2 s ) 
m 2 , < s < (1 - m 2 ) 


+0.3 
-0.3 
+0.3 
-0.3 


-6.08 
-6.08 
-0.391 
-0.391 


-24.0 
55.4 

0.276 
1.37 



Table 2: Values for the coefficients a(As)and /3(As). 
over a range (As): 

^l(As) = C 10 (a(As)C 9 + /3(As)) , (31) 

where 



a(As) 
/3(As) 

For a fixed value of A(As), one obtains hyperbolic curves in the C 9 -Ci 
plane; like the coefficients A, B and C, the parameters a and (3 are given in 
terms of the Wilson coefficients C\ ■ ■ ■ C§ and CV, and the kinematic region 
of s considered; their values are presented in Table |2|. 

Given the two experimental inputs, the branching fraction £>(As) and the 
corresponding asymmetry A(As), one obtains a fourth-order equation for the 
Wilson coefficients Cg and C\o, which admits in general four solutions. In 
figs. [l|-^ we plot the contours for a fixed value for the branching fraction 
B(As) and the FB asymmetry *4(As). Since £>(As) is an even function of 
Cio and ^4(As) is an odd one, we only plot positive values for C w . The 
asymmetry vanishes for C w = 0, but also for C 9 = — /3(As)/a(As). The 
two lines Cio = and Cg = — /5(As)/a(As) divide the C9-C10 plane into 
four quadrants, in which the asymmetry has a definite sign. Reflecting the 
hyperbolae on the line C\g = or Cg = —f3(As)/a(As) results in a sign 
change of the asymmetry. 

Figures [l] and |] show the contours in the Cg(fj,)-Cio(fj,) plane for the low- 
invariant-mass region 4m 2 e < s < m^, and figs. [3] and f| are for the high- 
invariant-mass region m^, < s < (1 — m s ) 2 . Figures |l| and ||| are obtained for 
CV(aO = 0-3? while figs. |] and [| are for C 7 (fi) = —0.3. The possible solutions 
for Cg and Cig are given by the intersections of the circle corresponding to 
the measured branching fraction and the hyperbola, corresponding to the 
measured asymmetry. Assuming the SM values for both £>(As) and ^4(As), 
one obtains the solid lines in figs. [L]-|]. The possible solutions in this case 
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8tt 2 f(m c /m b ) 

3a 2 1 
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dsu 2 {s)s (32) 
ds u 2 {s) \s Re Y(s) + 4C 7 (1 + m 2 )] (.33) 



are represented by solid dots (SM solutions) and solid squares (other non-SM 
possible solutions). 

From the figures one reads off that for the SM values of B and A one has 
more than one solution for the coefficients C 9 and C w , but the ambiguity 
may in general be resolved by measuring both the low and the high invariant 
mass regions. 

However, there is in principle also the possibility that the equations do 
not have a solution for Cg and C\q. This is the case, for example, when 
the asymmetry is large and the branching fraction small, in which case the 
hyperbola may not intersect with the corresponding circle any more. If this 
happens one has to conclude that the present analysis is not complete; in 
other words, the operator basis we started from is not complete and physics 
beyond the SM will be present in the form of additional operators such as 
right-handed currents. 

The upper and lower limits for the asymmetry as a function of the branch- 
ing fraction may be obtained analytically and are given by 



with 



A{As)± = Rcos9± (a(As)Rsin8± + y 



(34) 
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sin 9 
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a(As)R ^ \\ a 2 (As)R 2 
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+ 



(35) 

(36) 
(37) 



These boundaries are shown in fig. [5] for the low- and the high-invariant- 
mass region respectively. From these figures it is obvious that the bound is 
in fact non-trivial: it lies in any case far below the bound |.4| < B. 

If data become more precise one may in fact also expect a measurement 
of the spectrum and the asymmetry as a function of s. The spectrum itself 
is very sensitive to the values of the Wilson coefficients and to the sign of C-j. 
In fig. |6] we plot the various contributions to the spectrum, for positive and 
for negative Cj. 

In a similar way, it may become possible to measure also the differential 
asymmetry 



A{s) 



dz 



d 2 B 

dsdz 



dz 



d 2 B 

dsdz 



(38) 



which is also sensitive to the sign of C-j. The various contributions to A(s) 
are shown in fig. ^. 
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Figure 1: Contour plots of B(As) and A(As) in the C9-C10 plane for the low-invariant- 
mass region 4mj < s < m 2 j^ and C-j = 0.3. The circles correspond to fixed values of B: 
B = 5.6 x 10~ 6 (solid curve), B = 3.0 x 10~ 6 (long-dashed curve), B — 1.0 x 10~ 5 (short- 
dashed curve), B — 1.5 x 10~ 5 (dash-dotted curve). The left branches of the hyperbolae 
correspond to positive values of A: A = 1.7 x 10~ 7 (solid curve), A = 5.0 x 10 -8 (long- 
dashed curve), A = 5.0 x 10~ 7 (short-dashed curve), A = 1.0 x 10~ 6 (dash-dotted curve). 
The right branches of the hyperbolae correspond to negative values of A: A= — 1.41 ■ 10~ 8 
(solid curve), A = —5.0 ■ 10~ 9 (long-dashed curve), A = —3.0 ■ 10 -8 (short-dashed curve), 
A = —6.0 • 10 -8 (dash-dotted curve). For negative values of C 10 , the figure is simply 
reflected with A — > —A. The solid dot indicates the SM values for Cg and C\q. The solid 
square is another allowed solution resulting from the SM values of B and A. 
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Figure 2: Same as in fig. 1, but for C-j = —0.3. 

4 Model Predictions for the Wilson Coeffi- 
cients 

Once the parameters Cj, Cg, and Cio have been extracted from experimental 
data, a direct comparison with different theoretical predictions can be made. 
If a deviation from the SM result is observed, new models accounting for this 
discrepancy, can be looked for. If no such deviation is observed, the result 
can be used to set bounds on new physics. 

As already alluded to in the introduction, the operator basis considered 
here is not the most general one. For example, in the left-right symmetric 
(LR) models the basis (||) will have to be enlarged to incorporate the ex- 
tended SU(2)l x SU(2)u x U(l) gauge sector. The enlarged set contains 
the SM operators contained in T~i e ff above and another set of 10 operators 
obtained from this by flipping the chirality structure of the fermion fields 
Pl — > Pr, where Pl,r = (1± 7s)/2, increasing the number of independent 
Wilson coefficients to 20. 

Since we are restricting the operator basis to exclude ab initio LR-symmetric 
models, we give reasons why such extensions of the SM are unlikely to make 
a significant contribution to rare B decays. In the minimal LR-symmetric 
models, the right-handed CKM matrix is identical to, or the complex con- 
jugate of, the left-handed CKM matrix. However, in this case, there are 
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c 9 

Figure 3: Contour plots of B(As) and A(As) in the Cg-Cio plane for the high-invariant- 
mass region m^, < s < (1— m s ) 2 and for C 7 = 0.3. The circles correspond to fixed values of 
B: B = 2.56xl0~ 7 (solid curve), B = l.OxKT 7 (long-dashed curve), B = 5.0xl0~ 7 (short- 
dashed curve), B — 1.0 x 10~ 6 (dash-dotted curve). The left branches of the hyperbolae 
correspond to positive values of A: A — 1.41 x 10~ 8 (solid curve), A = 5.0 x 10~ 9 (long- 
dashed curve), A = 3.0 x 10 -8 (short-dashed curve), A = 6.0 x 10~ 8 (dash-dotted curve). 
The right branches of the hyperbolae correspond to negative values of A: A= — 1.41xl0~ 8 
(solid curve), A = —5. Ox 10~ 9 (long-dashed curve), A = — 3.0xl0 -8 (short-dashed curve), 
A = —6.0 x 10~ 8 (dash-dotted curve). For negative values of C w , the figure is simply 
reflected with A — > —A. The solid dot indicates the SM Values for Co, and Cio- The solid 
squares are other allowed solutions resulting from the SM values of B and A. 




Figure 4: Same as in fig. 3, but for CV = —0.3. 



rather strong lower limits on the mass of Wr, Mr > 1.5 — 2.5 TeV, arising 
from the condition that the short- distance contributions to the K^-Ks mass 



difference not exceed the experimental value [39|. In view of this restriction 
on Mr, it would be difficult to see tangible differences between the SM and 
the minimal LR models in rare B decays, much as has been argued for the 
mixing ratio x s and Xd j40] . These qualitative anticipations have been 
borne out by explicit calculations [|K], [Hj] . 



The situation with the non-minimal LR models is more involved [41]. 
While the bounds on ttir from ArriR can be evaded, the constraints from ex 



force rriR to be in excess of 30 TeV, in general [42]. So, we shall no longer en 



tertain the LR models, asserting that the present and impending constraints 
following from BiB — > X 8 j), Xd and Attir render the LR-symmetric-model 
effects in the FCNC semileptonic decays B — > X s £ + £~ insignificantly small 
and restrict the operator basis to the one given in (|4]). In this section, we con- 
sider as illustrative examples two popular extensions of the SM (the MSSM 
and the 2HDM) and study how the predictions for the Wilson coefficients 
C 7 , Cg, and C w are altered in these models, compared to the SM. 

Supersymmetric models have a new source of flavour-changing neutral 
currents because of the quark-squark-gluino vertex, which is, in general, not 
diagonal in generation space [131. If the generation mixings of such ver- 



tices were arbitrarily large, they would lead to phenomenologically unac- 




Figure 5: The bounds for the partially integrated asymmetry 
A(As) as a function of the partial dilepton branching ratio 
B(As). The upper plot is for the "low-invariant-mass region" 
and the lower one is for the "high-invariant-mass region" as 
defined in the text. The boundary demarcated by the solid 
curves are for the positive value C 7 = 0.3, the long-dashed 
curves correspond to the negative value CV = —0.3. The 
dotted line is the trivial bound |*4| < B. The solid dot denotes 
the SM value. 




Figure 6: The dependence of the invariant-mass spectrum 
on the Wilson coefficients. Solid line: SM. Long-dashed line: 
Cj — > — Cj, with other coefficients retaining their SM values. 
Short-dashed line: The contribution of C w only. Dotted line: 
Cio = 0, with other coefficients retaining their SM values. 
Dash-dotted line: same as for the dotted one, but with C 7 = 
—0.3. The vertical lines indicate the location of the J/\I/ and 
\P' resonances. 
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Figure 7: The dependence of the differential FB asymmetry 
on the Wilson coefficients. Solid line: SM. Long-dashed line: 
C-j — > — CV, with the other parameters retaining their SM 
values. The vertical lines indicate the location of the J/^ 
and ^' resonances. 
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Figure 8: The dependence of CV on the charged Higgs mass 
for V2/V1 = 1 (solid line), Vijv\ — 1 (dashed line), v<ijv\ = 10 
(dot-dashed line). The band between the two curves reflects 
the QCD uncertainty (see text). 



ceptable flavour violation in the K°-K° system. However, in the minimal 
super symmetric model, flavour universality of the supersymmetry-breaking 
terms is usually assumed to hold at the grand-unification scale. In this 
case, a non- vanishing (and calculable) generation mixing in the quark-squark- 
gluino vertex is induced only by the effect of renormalization from the grand- 
unification scale to the weak scale. Given the present experimental limit on 
the gluino mass, this effect leads only to very small contributions to flavour- 
changing ^decays [0 and it will be neglected in our analysis. Gluino-mediated 
flavour- changing neutral currents may however play an important role in non- 
minimal supersymmetric models (see [44]] and references therein). 



Once the gluino contributions (as well as the analogous ones from neu- 
tralino exchange) are neglected, the flavour violation in the model is com- 
pletely specified by the familiar CKM matrix. The one-loop supersymmetric 
corrections to the Wilson coefficients C 7 , C 9 , and C w are given by two classes 
of diagrams: charged-Higgs exchange and chargino exchange. Their analyti- 
cal expressions can be found in ref. [J/J, and will not be repeated here. 

The charged-Higgs contribution is specified by two input parameters: the 
charged-Higgs mass (m H +) and the ratio of Higgs vacuum expectation values 
{v-ijv\ = tan/3). This contribution is interesting in itself, since it corresponds 
to a well-defined SM extension: the 2DHM. The interaction between the 
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Figure 9: The dependence of Cg on the charged Higgs mass. 
The notation is the same as in fig. || 



charged Higgs and quarks is given by: 

C = ^£^ H+U [A U M U V^—^- + A d VM d ^^j d + h.c, (39) 

where M U}d are the up and down quark mass matrices and V is the CKM 
matrix. If both up and down quarks get masses from the same Higgs doublet 
(this case is usually referred to as Model I), then 

A u = -A d = l/tan(3. (40) 

In the supersymmetric case, two different Higgs doublets couple separately 
to up and down quarks (Model II) and 

A u = 1/A d = 1/ tan/3. (41) 

Let us first consider the case relevant to supersymmetry, i.e. Model II. 
Figures §, § and [LC] show the charged-Higgs effects on CV, Cg, and Cio for 
m t = 174 GeV, at leading order in QCD. For each value of v 2 /vi, two curves 
for r] = a s (M z )/a s (n) equal to 0.68 and 0.41 are shown. The band within 
the two lines reflects the QCD uncertainties, since the two chosen values of 
rj correspond to the uncertainties in a s {Mz) (= 0.118 ± 0.006) |46| and fi 
(nib/2 < fi < 2mb)- Recently the known part of the next-to- leading correc- 
tions has been included in the analysis of b —>■ s'j ^5|. They have the effect 
of reducing the scale dependence, and of decreasing considerably the value of 
Cj with respect to the leading-order calculation shown in fig. ||. We expect 
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Figure 10: The dependence of C\g on the charged Higgs 
mass. The notation is the same as in fig. [| 



that we can safely ignore the next-to-leading corrections to Cg, because the 
QCD running is much less important here than in Cj. Note that Cio does 
not get renormalized under QCD. 

The constraints imposed by supersymmetry on the scalar potential re- 
quire tan/? > 1 and m H + > my/- Notice that, as tan/? becomes large, Cg 
and Cio rapidly converge to the SM value (as tan -2 /?), while the charged- 
Higgs contribution to Cj remains non- vanishing even in the limit tan j3 — > oo. 
In the case of a non-supersymmetric Model II, the experimental limit on 
BR{b — > 57) rules out a large region of the tan/3-m^ parameter space j|, §]. 
As a consequence, no significant deviation from the SM values of Cg and 
Cio can be expected. In supersymmetry, however, the chargino can largely 
compensate for the charged-Higgs contribution to Cj [§], and more decisive 
differences from the SM are foreseeable, as discussed below. 

In the case of Model I, the expressions for Cg and Cio are still the same as 
in Model II. However C7 is modified and the constraint from BR(b — > sj) on 
the parameter space is much weaker. For tan /3 < 1, C7 can become negative. 
For small mjj+ and tan (3, it is possible to reach values for which Cj equals, 
in absolute value, the SM prediction, but is opposite in sign. The region 
of parameter space where this happens is barely allowed by the constraints 
from 6k and B°-B° mixing | 47f , but leads to unacceptable corrections to the 
Z° — > bb width ESI and is thus ruled out experimentally. 



In addition to the diagrams with charged-Higgs exchange, the MSSM 
leads also to chargino-mediated diagrams. The chargino contribution is spec- 



ified by six parameters. Three of them enter the 2x2 chargino mass matrix: 



m x+ = ( % . ™wy/28m{l\ (42) 
x \mw\ 2 cos p jj J 

Following standard notations, we call tan j3 the ratio of vacuum expectation 
values, the same that appears also in the charged-Higgs sector, and M, \x the 
gaugino and higgsino mass parameters, subject to the constraint that the 
lightest chargino mass satisfies the LEP bound, m+ > 45 GeV. The squark 
masses 

m?2 = m 2 + m 2 q ± Amm q (43) 

contain two additional free parameters besides the known mass of the cor- 
responding quark m q : a common supersymmetry-breaking mass m and the 
coefficient A. The parameter A contains all the information (both from the 
/i-term and the trilinear term) of left-right squark mixing and it is constrained 
by the requirement that the lightest stop is not produced at LEP, m 4 ~ > 45 
GeV []. We also take into account the CDF limit on squark masses and 



impose m > 126 GeV. In eq. fl43|) we have made the simplifying assumption 
that the supersymmetry-breaking left- and right-squark masses are equal. 
The last parameter included in our analysis is a common mass mj for slep- 
tons, all taken to be degenerate in mass, with the constraint mj > 45 GeV. 
Therefore the version of the MSSM we are considering is defined in terms of 
seven free parameters. 

We have computed the Wilson coefficients in the MSSM and then varied 
the seven above-defined parameters in the experimentally allowed region. 
The results of our analysis are presented in fig. OL which shows the regions 
of the C9-C10 plane allowed by possible choices of the MSSM parameters. 
The upper plot of fig. |TT| corresponds to parameters which give rise to positive 
(same sign as in the SM) values of C7, consistent with experimental results on 
b — > 57 (0.19 < C-j < 0.32), while the lower plot corresponds to values of C7 
with opposite sign (—0.32 < C 7 < —0.19). We also show how our results are 
affected by an improvement in the experimental limits on supersymmetric 



particle masses, as can be expected from the Tevatron and LEP 200. Fig. [11 
also shows the Cg— C\q regions allowed by the MSSM if the further constraints 
rriH+ > 150 GeV, m ( -, m x +, mj > 100 GeV are imposed. 



The regions shown in fig. |TT| illustrate the typical trend of the super- 
symmetric corrections. If supersymmetric particles exist at low energies, we 
can expect larger values of C w and smaller (negative) values of C 9 than 
those predicted by the SM. This is the general feature, although the exact 
boundaries of the allowed regions depend on the particular model-dependent 
assumptions one prefers to use. For instance, the allowed regions can be 
slightly expanded if one introduces more general stop mixings than those 

4 For a particular choice of the mixing between the two stop states, the light stop can 
be decoupled from the Z° and the LEP bound would not apply. We disregard here this 
possibility. 
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Figure 11: The region in the Cg-C w plane obtained by 
varying the MSSM parameters. The upper (lower) plot cor- 
responds to solutions that satisfy the b —>■ 57 experimental 
constraint with positive (negative) Cj given in eq.(|) and the 
present bounds (m#+ > 80 GeV, m t , m x +, rfi£ > 45 GeV). The 
smaller areas limited by the short-dashed line correspond to 
the region of the MSSM parameter space that will survive an 
unsuccessful search for supersymmetry at the Tevatron and 
LEP 200 {m H+ > 150 GeV, m h m x+ , m t > 100 GeV). 



considered here. On the other hand, if further constraints are imposed on 
the model (such as particular boundary conditions at the GUT scale, radia- 
tive symmetry-breaking, etc.), some of the parameters which in our analysis 
were taken as independent variables become related with each other, and the 
allowed region may somewhat shrink. However, the most interesting feature 
of supersymmetry is that solutions with negative values of C-j are possible 
and are still consistent with present data. Moreover, values of the other two 
coefficients Cg and C\q sufficiently different from the SM are allowed, leading 
to measurable differences in the decay rates and distributions of B — > X s £ + £~ 
and B s -> £+£-. 



5 Concluding Remarks 

The branching ratios B(B — > K*j) and B — > X s 7 have provided first 
measurements of the strength of the effective magnetic moment operator 
rribSLO nvbnF^ '. The inferred value of the Wilson coefficient |CV(m{,)| =0.19 
- 0.32 from present data is in broad agreement with the SM estimate of the 
same. This, in turn, has led to rather stringent bounds on the parameters of a 
number of models. The comparison between the SM and experiment has be- 
come more precise with the measurement of the inclusive rate B(B — > X s ^f). 
There is an overriding need to firm up theoretical estimates in the SM by 
calculating the next-to-leading-order effects. 

While the first step in the measurements of rare B decays has been made, 
we expect this field to undergo a qualitative change as the anticipated exper- 
imental facilities for B physics take shape and start producing B hadrons in 
large quantities. We have in mind here a steady consolidation of the CLEO 
and CDF data, with big strides expected to be made at the threshold B 
factories at SLAC and KEK and experimental facilities with proton beams 
such as HERA-B and in particular the LHC. Among others, we expect mea- 
surements of the CKM-suppressed radiative decays B — > p7 and B — > uj, 
and the FCNC semileptonic and leptonic decays B — > X s £ + £~ (and the 
corresponding exclusive decays B — > (K, K*)£ + £~) and B° s — > £ + £~ at these 
facilities. Also, the nature of rare 5-decay measurements will evolve in time 
from being exploratory to becoming rather precise. This has encouraged us 
to propose undertaking a more ambitious programme of extracting from data 
model-independent quantities - very much along the same lines as was carried 
out for the precision electroweak analysis of the LEP, SLC and low-energy 
data. 

The electroweak precision tests, at LEP and elsewhere, have concentrated 
on the self-energies of the electroweak gauge bosons (7, Z° and W ), with 
possible deviations from the SM expressed in terms of a limited number of 



parameters [50[ 51]. The main interest in rare B decays is to measure the 



effective FCNC vertices to test the SM precisely and search for new physics. 
We have argued how to parametrize these vertices through a limited num- 
ber of effective parameters, which govern the rates and shapes (differential 
distributions) in rare B decays B — > X s 7, B — > X s £ + £~and B s — > £ + £~ . 
The search for physics beyond the SM in these decays can be carried out 
in terms of three effective parameters, which can then be interpreted in a 
large class of models. The presence of non-SM physics may manifest itself 
by distorting the differential distributions in B — > X s £ + £~ . Some possible 
examples of such distortions have been worked out. While the analysis pre- 
sented here covers a large class of models, we have also presented profiles of 
the Wilson coefficients in the best-motivated extensions of the SM, namely 
the MSSM. We have pointed out that there exist two distinct solutions cor- 
responding to the negative and positive values for CV(/i), which cannot be 
distinguished from the data on B(B — > Xg'y), but they give rise to very dif- 
ferent distributions in the decay B — > X s £ + £~. Moreover, the coefficients of 
the operators Cg(/x) and Cio(n) f° r the MSSM models may attain values suf- 
ficiently different from the corresponding SM values. An encouraging result 
of our analysis is that rare B decays B — > X s ^y and B — > X s £ + £~ have a 
discovery potential much beyond direct searches at these facilities. Finally, 
we remark that the purely leptonic decay modes B® — > £ + £~ can also be 
used to further pin down the coefficient of the operator Oiq involving the 
axial- vector leptonic current, i.e. |Cio(/i)|, although this may require much 
larger _B-hadron statistics. Precision measurements in rare B decays must 
be pursued vigorously. 
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Appendix: The effective Hamiltonian for B — > 

X s7 and B -> X s tt 

The effective Hamiltonian is 

AC 10 

H eff = i^-KVtb^Cii^Oiifi), (44) 

V2 i=1 

where the operator basis is chosen to be 

Ox = (s La 7 /i fe La )(c L/3 7 At c L/3 ) (45) 

02 = (sLa7nbLp)(cLp7»CLa) (46) 
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2 1 111 2 
The Wilson coefficients of the operators are given by the renormalization 
group evolution 
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To leading logarithmic level, the QCD beta function (3(g) is given by 

Pig) = -^tJ^ with Po = n- \f 



(60) 



and 'j(g) is the anomalous dimension matrix, which is, to leading logarithmic 
acuracy, given by 



lis) = -7o 
Here 70 is a 10 x 10 matrix given by 

7o = 



9' 



16tt 2 



(61) 



9 

— z 


D 


n 
u 


u 


n 

u 


n 
u 


n 

u 


q 
o 


16 

3 






u 





— z 


2 
9 


2 
3 


2 
9 


2 
3 


416 
81 


70 
27 


16 
9 






n 
u 








22 
9 


22 
3 


4 
9 


4 
3 


464 
81 


^+3/ 


-T?(«-3 















6-1/ 


-2 + !/ 


~-f 


-f 

3-' 


416 232 d 
81 U 81 " 


6+1/ 


-f (--1 




3) 

















2 


-6 


32 


-T-3/ 




d 
2 


) 











--f 


-f 


9 J 


-16+|/ 




-4-W/ 




d 
2 


) 























32 
3 

































32 
9 


28 
3 

















































































Using as initial condition Cj(Mw) = for j = 1, 3, • • • 6 we obtain, for the 
solution of the renormalization group flow 

CM = \C 2 (M W ) ( V ^ 3 - V -^) (62) 

CM = \c 2 {M w )(^ + r 12/2 *) (63) 

ri(..\ - r< ( M \ ( n m 1 0^0.8994 , 1 -12/23 n 1 /ino^-0.4230 , n nnK/1^0.1456 /^a 



CM = C , 2 (M H ,)(-0.0112r ? a8994 + ^- 12 / 23 - 0.1403r ? - a4230 + 0.005V' 1456 (64) 

- 0.07147? 6 / 23 + 0.0509?? - 4086 ) 
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and we use C 2 (Mw) = — 1. 

In terms of this effective Hamiltonian, the amplitude for b — > s£ + £~ be- 
comes 
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where = Pi + p% is the momentum transferred to the leptons. The ef- 
fective coefficient Cg e // contains the contributions from the one-loop matrix 
elements of 0\,...,Oq and is given by (||). This definition involves the one- 
loop function g, which is given by 
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